N86 - 30185 


MAGNETIC BEARINGS FOR VIBRATION CONTROL 
G. Schweitzer 

Swiss Federal Institute of Technology (ETH) 
Zurich, Switzerland 


The paper presents a survey on the research of the Institute of Mechanics of the 
ETH in the field of vibration control with magnetic bearings. It shows a method for 
modelling an elastic rotor so that it can be represented by a low order model amen- 
able to control techniques, it deals with the control law and spill-over effects, and 
it discusses experimental results for an active resonance damper. 


INTRODUCTION 


Magnetic bearings can exert forces on a body without any physical contact. This 
makes it a very useful device to influence the position and the motion especially of 
a spinning rotor. Magnetic bearings are quite capable of supporting even a heavy ro- 
tor. Its application to turbomachinery, machine tools, and in the vacuum techniques 
is described for example by Habermann (ref. l). Such a bearing, or more precisely 
such a bearing system, can be designed in a systematic way (ref. 2). By suitably de- 
signing the control loop the magnetic actuator can be adjusted to a variety of appli- 
cations and dynamic requirements. Primary design goals have been to support a rigid 
rotor, because the rigidity of the rotor facilitates the control design essentially. 
The magnetic forces can be made to be a function of the rotor motion in such a way 
that the actuator usually has spring and damper characteristics which suitably depend 
on the excitation frequency. 

The freedom in assigning dynamic characteristics to the actuator can be used not 
only for supporting a rigid rotor but for vibration control as well. Some efforts 
have already been undertaken to control various kinds of vibrations. Pietruszka and 
Wagner (ref. 3) show how an unbalanced rigid rotor can be made to rotate about its 
principal, axis of inertia, thus avoiding vibrational unbalance forces on the magnetic 
bearings. In reference L magnetic bearings are used to shield a measuring platform 
from residual vibrations of its base. Gondhalekar et al. (ref. 5,6) investigate vi- 
bration control problems of an elastic rotor. Even the active damping of self-excited 
vibrations by internal friction and of parametric vibrations due to rotor asymmetry 
have been looked into. It shows that pragmatic solutions for individual problems can 
be found, general methods and the answers to some basic questions, however, are not 
yet available. 

This paper will give a survey on our research in the field of vibration control. 
At first it presents the modelling of an elastic rotor so that the resulting low 
order model is amenable to control techniques . Another aspect concerns the control 
law and the effects of spill-over, and finally design requirements and experimental 
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results for an active ''resonance damper are discussed. The survey is based on the work 
of the research group of the Institute of Mechanics at the ETH, especially on the 
contributions of Bucher, Salm and Traxler. 


FUNCTIONAL PRINCIPLE 



Fig. 1: Magnetic Suspension 



At first the functional principle 
of the magnetic suspension will be ex- 
plained shortly and demonstrated by 
examples . 

Figure 1 shows a hovering shaft S. 
Its position is measured by a sensor, 
the measured signal is processed in a 
controller, the controller output con- 
trols via an amplifier the current 
through the coils of the magnet M in 
such a way that the magnetic force 
stabilizes and firmly keeps the shaft 
in its hovering state. This simple 
example shows that obviously the size 
of the magnetic force depends on the 
size of the air gap, and that the dy- 
namics of the suspension, its stiffness 
and damping, to a large extent is de- 
termined by the design of the control- 
ler. 


For suspending a full rotor the 
simple loop will not be Sufficient. 
Figure 2 shows the radial suspension 
of a rigid rotor. For each degree of 
freedom a magnetic actuator has to be 
controlled independently. In general 
each bearing force will depend on each 
sensor signal, leading to a typical 
multivariable control. 


Fig. 2: 


Block-diagram for the radial 
suspension of a rotor in one plane 
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Figure 3 shows the hardware setup 
for such a bearing system. The rotor 
has a length of about 1 m, a diameter 
of 12 cm and a mass of 12 kg. The air 
gap of 1 cm is extremely large for 
technical purposes. The device was 
used for an exposition. The signals 
of the optical CCD-sensors are pro- 
cessed by a microprocessor and fed to 
switched power amplifiers. 


Fig. 3: Magnetic bearing system, front 

view with control unit , power 
supply and drive 


MODELLING OF AN ELASTIC ROTOR 


For most of the current applications the rotor is considered to be rigid. The 
derivation of its equations of motion under the action of magnetic control forces 
shouldn't cause any difficulties. For an elastic structure the modelling procedure 
is somewhat more complex, as will be outlined in the following. 

The model of the elastic rotor has to be such that it can be incorporated into 
the model of the closed control loop, and it has to be of low order so that the com- 
plexity of designing the control law will be reduced and a simulation of the closed 
loop becomes reasonable. 

For the derivation of such a model, following Bucher (ref. 7), let us start 
with a finite element model of the elastic rotor 

(!) ^ £ + A x 4 + ^ £. * R 

where q is a n x 1 vector of generalized displacements, £ is the vector of generalized 
forces, and the A£ are structural coefficient matrices, characterizing inertia and 
elasticity as well as gyroscopic and nonconservative properties. Equation 1 may be 
transformed into the complex frequency domain and expressed by 

(2) [s 2 Ag + s A x + Aj Q(s) = P(s) 

The matrix in brackets is termed dynamical stiffness matrix and its inverse is the 
dynamical flexibility matrix 
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(3) 


H(s) = [s 2 A 2 + sA^ + aJ _1 

The elements of H(s), the dynamical flexibility transfer functions, are well-known in 
measurement and modal analysis techniques. A general way of reducing the large set of 
finite element equations is to truncate the modal representation of equation 1. In 
order to find that representation we have to solve the eigenvalue problem of equation 
1. For the sake of simplicity - all the subsequent derivations can be extended to the 
general system (equ. 1) as well (ref. 7) - let us assume a simple elastic structure 

(U) M q + K q = £ 

From the solution of the eigenvalue problem we obtain a set of eigenvalues uj^,...u) n 
and the corresponding real normalized eigenvectors uj^ with the modal matrix 
IJ = ( U1 , . • • »Ujj » • • • »Ujl ) . 

Then the modal expansion of q 

(5) q * I u^. z k a U. z , k = 1, . . . ,n 

leads to the modal representation of equation k 

(6) Mz_+Kz_=U^p 

■with M « U T U = diag(m^), K = U T U = diag (k^) , a) 2 = k^/m^ 

This modal representation is a time domain representation of the partial fraction ex- 
pansion of the corresponding flexibility matrix H(s). Translating equations 5 and 6 
into the s-domain, yields 

T 

2 -IT n u* 

(7) I(s) - U(M s 2 + K) 1 U T = I — — 

k=l m^( s +(jJ j s ) 

For the design of the controller , . a model of the flexible structure, that relates 
the displacements and velocities at the location of the sensors to the excitation for- 
ces of the actuators, is required. By analyzing the closed loop system with respect to 
the actuator and sensor coupling points, it can be checked whether the design objec- 
tives have been reached, i.e., whether the system is stable and shows up with good 
damping performance. Also, test measurements for the closed loop system are easily 
carried out by exciting the structure with the available actuators. Hence, for both 
control design and closed loop analysis , we may confine the selection of coordinates 
to the coupling points of the actuators and sensors with the structure . A modal re- 
presentation of the dynamical flexibility matrix of the elastic structure, with re- 
spect to the cited coordinates, provides the desired description, as will be shown 
subsequently. 

The n^ actuator forces form the vector w, the n^ measured displacements are given 
by £ , and the n T measured velocities by y v . These vectors are related by appropriate 
incidence matrices to the nodal forces and displacements by 
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( 8 ) 


T 

^41 


*D = — D 


y = T 

v "v — 


T , T T, 
<1 L = %, 


The elements of the incidence matrices are zero or one, when discrete actuators and 
sensors are located at nodal points. We assume that we can truncate the number of 
modes retaining only m essential modes, so that instead of equation 5 we now have 


(9) q “ £ u . z . = U z , i = 1 , . . . ,m , m<n 

. — l l - m - m 

l 

with the truncated n x m modal matrix U and the truncated modal vector z . This leads 

. — m — m 

to the truncated set of modal equations 


( 10 ) 


M z + K z 
— m - m - m ~ m 


T 

= U 
— m 


P 


With the subset of retained coordinates (equ. 8) we have in the s-domain a transfer 
matrix representation 


( 11 ) 


Y( s ) = W( s ) Hj, 


U (s 2 M + K ) 1 U T T A 
— D ~m -m — m ~rn —A 


s T U (s 2 M + K ) 1 U T T a 
— v -m -m -m -m — A 


Here the transfer matrix Hq is a generalized dynamical flexibility matrix, since it 
relates a force excitation (actuator forces) to nodal displacements and velocities 
( sensor signals ) . 


By introducing the state vector 




,T 

z 


— m 


) 


an equivalent state space representation is obtained 


(12) 


A = 


x = A x + B w 
X = £ x 


0 

-m' 1 K 

— m — m 


E 

0 


B = 

0 

C = 


T 

U T a 
— m —A 



T- u 


0 


0 

T U 

— y -m 


All of the representations (equ. 10,11,12) will be suitable for the control design. 
The method can be extended to include other structural elements as well, i.e. found- 
ation dynamics, by applying the building block approach. Problems left over and to be 
discussed in the next section are that of truncation and the control, design. 


CONTROL OF AN ELASTIC ROTOR 

When the control design is based on a reduced-order model of the rotor (equ. 10 ) 
the question arises whether the real rotor will indeed be controlled correctly by 
such a simplified controller. Actually there may be detrimental effects, and we can 
explain and classify them in the following way (ref. 6). 
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Let us partition the high-dimensional modal coordinate vector _z into two parts. 
The important part with the low dimension n m will he the one to he controlled, the 
other residual part of dimension n r is the one to he neglected. Then equation 6 and 
the measurement equations 8 can he arranged in the following form 


- 
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0 


z 


K 

0 


z 


u T 

-m 



— m 
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. — r 


= T_ U z + T„ U 
— D - m ~ m — D — r 


*v 


■ T U 
—v — m 
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+ T 
-v 


U 


For design purposes it will he assumed that the controller is based on the re- 
duced-order model only, being a subset of equations 13 


(14) 


M z + K z a U T p = U T T a w 

— m — m — m — m — m z. — m —A ~~ 


Zjte - 1 D u, ^ . 


*vm 


T U z 

— v -m ~m 


The relation between the control vector w and the measurement follows from the chosen 
control law, which preferably will be a linear one and can be determined in the usual 
ways. Such a control design will result in the desired and "good" control for the 


reduced-order model (equ.l4), but when it is 


Rduator Rotor Senior 



Fig. 4: Reduced order control applied to 

the real full-order system, de- 
monstrating the observation spill- 
over (OS) and the control spill- 
over (CS) 


applied to the real full-order-system 
(equ. 13) the system qualities can 
alter essentially (fig. 4). In reality 
the measurements £ do not only consist 
of the modelled part but they also 
depend on the residual vector z , 
causing the socalled "observation 
spillover" . 

Furthermore the control vector w 
does not only act on the modelled part 
but on the real full-order system it- 
self which obviously contains the un- 
modelled part as well. This influence 
is called "control spillover". These 
spillover terms can change and deteri- 
orate the behaviour of the real system 
and even destabilize it. 

The objective of the control de- 
sign now is to derive the control on 
the basis of the reduced-order system 
without exact knowledge about the re- 
sidual, unmodelled part, and to pro- 
duce a design that provides suspension 
and damping for the mechanical system 
and that is easy to implement as well. 
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A direct output feedback is known to have these properties and to show robust- 
ness qualities with respect to parameter errors and truncation effects (ref. 8). Then 
the control forces depend linearly on the displacement and velocity measurements 


(15) 


2- -2d*d '2^ 


with the yet unknown gain matrices G- . Introducing this control law and the measure- 
ment equations 8 into the full-order system equations 6 leads to the full-order 
closed-loop system (FOCL) 

(16) M ’z + U T T„ G T U z + (K + U T T A U)z = 0 

In order to avoid the spillover effects the solutions of equation l6 should be stable 
Then the FOCL is stable despite the fact that its gain matrices will be based on the 
reduced-order model. The stability conditions, derived from definiteness properties 
of the coefficient matrices (ref. 6), are certainly fullfilled, 

. T 

- if the sensors and actuators are arranged in collocated pairs, i.e. if T. = T^ = T_^ 

And when the actuators have to support the rotor additionally to controlling its 
vibrations the number of reasonably placed sensors and actuators has to be at least 
equal to the number of not stable rigid body modes. 

- if the gain matrices G. are symmetric and positive definite. 

The procedure how to construct such gain matrices will be explained shortly. The re- 
duced-order model, including the rigid-body modes, and all other vibration modes that 
should be controlled, is given by equation lU . Introducing the control equation 


— — D ^-Dm ^vm 

into equation lU leads to the closed-loop equation of the reduced-order model 


(IT) 

( 18 ) 


M z+D z + (K+S)z = 0 

-m — -m — -m -m — 


D = U T. G T U , 
— m — m —A ~v ~ v — m 


S = U T G T U 
— m - m —A — D — D — m 


The terms D , are the additional modal damping and stiffness introduced by the con 
trol feedback. The required symmetry and definiteness properties of the yet unknown 
gains G^ will be guaranteed by a design approach of Salm (ref. 6). First we assign 
desired additional modad damping and stiffness values and thereby choose the diagonal 
matrices 1^ and Sj^, and then we solve equations 18 for the unknown gains. However, 
as most often the number of n^, n-p of actuators and sensors will be smaller than the. 
number m of modelled modes, it is not possible to solve equation 18 uniquely. An 
approximation for the gains G_^, in the sense of minimized quadratic differences using 
Pseudo-Inverses, is feasible and leads to the following design rules: 

The reduced-order control of large, and even infinite dimensional vibration 
systems will result in a robust and certainly stable FOCL, if 
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- the vibration system is not unstable (rigid body modes are allowed), 

- the reduced-order model contains at least all rigid-body modes , 

- the sensors and actuators are collocated and their number is at least equal to the 
number of rigid-body modes. 

This control approach includes the technically most interesting case where a real 
flexible rotor is suspended by actuators which at the same time have to control its 
elastic vibrations. An example and experimental results will be shown in the next 
section. 


EXAMPLES AND EXPERIMENTAL RESULTS 


Active suspension and vibration control of a flexible beam (ref. 6). The trans- 
verse vibrations of a flexible beam, supported and controlled by actuators on both 
ends, are described by the modal equation 6. The mode shapes (equ. 5) used in this 
case are the free-free modes, the rigid body modes are included in the modal trans- 
formation matrix U. 


As a design goal the two actuators will 
and additionally two elastic modes should be 



have to control the two rigid-body modes , 
strongly damped, too. Therefore the re- 
duced-order model is chosen to con- 
sist of the two rigid-body modes and 
the first two elastic modes. Dis- 
placement and velocity signals are 
coming from sensors at the actuator 
locations. The gain matrices for 
the direct -output control are derived 
from equation 18, where in our case 
the assigned modal stiffness and 
damping, that we want to obtain, is 
characterized by 

2 

S = Diag w . D = Diag (2D w . ) 

— m mi — in mi 

The stiffness of the suspension was 
assigned by prescribing 0 )q, u)2* The 
damping value D was varied between 
0 and 1. Several general and anti- 
cipated results could be corroborated 
numerically: 


Fig. 5: Typical eigenvalue-curves of the 

controlled full-order system with 
increasing damping values D 
assigned to the modes 


- For"low" bearing stiffness (Fig. 5a) 
the eigenvalues of the controlled 
system have a typical behaviour 
when we increase the assigned damp- 
ing value D. For very high damping 


the modal frequencies of a fixed-fixed beam are approximated. 


- For a "stiff” bearing design (Fig. 5b), where we want the bearing stiffness to be 
near the first elastic frequency of the beam, the behaviour is similar, but some 
eigenvalues now go upwards to higher eigenfrequencies of the fixed-fixed beam. The 
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poor controllability of these modes comes from the fact that high control forces 
tend to "clamp" the beam, and to shift modal nodes- into the actuator location it- 
self. 

- Limiting values for the attainable stiffness are given by the lowest frequencies of 
an equivalent beam with fixed ends. 

- The control design always produces a stable system. It provides stability and dam- 
ping for the modelled and the unmodelled modes. 



Fig. 6: Rotor model, demonstrating the 

parameter deviations £,? and 
61 from a reference position 



Fig. 7: Experimental setup 


Vibration damping and spill-over effects . 

A simple rotor model (Fig. 6) demonstrates 
the active damping of resonance vibrations 
(ref. 6). The experimental setup (fig. 7) 
consists of a conventionally supported, 
flexible rotor shaft with a disc in the 
"middle" of the shaft. The "exact" model for 
the rotor of figure 6 considers the imper- 
fection parameters £,£,61. Now the control 
design is based on a reduced-order model, 
taking only the displacements Uy, u z into 
account. The motion of the "exact" rotor, 
however, will be affected by the coupling 
between translation and inclination of the 
disc. It is the imperfection parameter £ that 
causes observation spillover, and £ pro- 
vides control spillover. The eigenvalue 
shifts caused by these spillover effects are 
demonstrated in figure 8. The unmodelled 
modes can even become unstable. When actua- 
tor and sensor are collocated (e = £) then 
the "exact" model is stable even with the 
"simple" control. The calculations have 
been confirmed by experiments. Near the 
stability boundary however the above model- 
ling was not sufficient: the small imper- 
fections due to the phase/frequency response 
of the sensor had to be taken into account 
(ref. 7). 

When the elastic rotor and the control 
are modelled according to the above rules , 
the active damper can reach good performance 
Measurement results (fig. 9) demonstrate 
the considerable reduction of the resonance 
amplitudes, for the experimental setup. 
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CONCLUSION 



Fig. 8: Eigenvalue shift of the un- 
modelled modes depending on the 
parameters e and B,. Without 
control the eigenvalue is at A. 



O 40 

Fr«qu«nz HZ 


In order to control vibrations in flex- 
ible mechanical systems and to solve special 
suspension tasks it can be advantageous to 
use magnetic actuators. For the design of the 
multivariable control loop a suitable model 
of the mechanical system is derived. Starting 
from a FE-model of high order the reduced 
order model finally only retains the actua- 
tor and sensor coordinates and the modes 
considered to be important. In the control 
design the spillover aspects are discussed. 

A method for determining a control based on 
the reduced-order model and robust with re- 
spect to the number of unmodelled modes is 
demonstrated. Examples show the active sus- 
pension of a flexible beam, spillover effects, 
and measured resonance curves for an active- 
ly damped rotor running through its critical 
speed. 


Fig. 9: Measured amplitude/ frequency 
response, (a) uncontrolled, 
(b) controlled 
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